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I present a partly pedagogic discussion of the Gel'fand- Yaglom for- 
mula for the functional determinant of a linear, one-dimensional, sec- 
ond order difference operator, in the simplest settings. The formula 
is a textbook one in discrete Sturm-Liouville theory and orthogonal 
polynomials. A two by two matrix approach is developed and applied 
to Robin boundary conditions. Euler-Rayleigh sums of eigenvalues 
are computed. A delta potential is introduced as a simple, non-trivial 
example and extended, in an appendix, to the general case. The con- 
tinuum limit is considered in a non-rigorous way and a rough com- 
parison with zeta regularised values is made. Vacuum energies are 
also considered in the free case. Chebyshev polynomials act as free 
propagators and their properties are developed using the two-matrix 
formulation, which appears to be novel and has some advantages. A 
trace formula, rather than the more usual determinant one, is derived 
for the Gel'fand- Yaglom function. 
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1. Introduction 

Finite structures are very common in science either as approximations to some 
continuous arrangement, perhaps for numerical purposes, or because of some in- 
herent discreteness or, again, for regularisation. They have also gained a certain 
currency in elementary particle models. 

In this communication I wish to make some rather elementary computations of 
one or two quantum field theory quantities using finite difference notions. 2 I restrict 
myself to the simplest one-dimensional systems i.e. fields on the interval or circle. 

Although these have been discussed, almost ad nauseam, I could not find this 
particular development completely in the literature. The interesting work by Actor, 
Bender and Reingruber, [4] , contains a detailed treatment of the Casimir effect on 
the lattice, and, while I cannot add too much to their extensive results, I will 
recover some of the formulae for completeness. I will also compute the discrete 
determinants for the free field case and I will include a mass here too. Although the 
explicit results are rather trivial, and just examples of general expressions, I believe 
they have some didactic merit. As something more substantial, I also treat the case 
of a delta potential. 

Functional determinants appear in many areas and their computation is im- 
portant physically. An early method is the Gel'fand-Yaglom technique which is a 
means of finding an operator determinant without knowing the eigenvalues explic- 
itly. The continuum case (originating with Gel'fand and Yaglom, [5], and Levit 
and Smilansky, [6]) has been analysed fairly extensively. The work by Kirsten and 
McKane, [7] , contains a brief historical survey plus a contour integral proof of the 
theorem and a discussion of the zero mode problem. In the quantum field theory 
context, Dunne, [8], can be consulted for orientation and further references. 

There has been less work on discrete systems, although there is a considerable 
body of work concerned with graphs, which, though relevant, I will not consider, 
per se. Very general theorems have been derived by Forman, [9], for the situation 
when a potential is present. He proves and employs a discrete Gel'fand-Yaglom 
theorem. In the following sections, I give a simple justification of the formulae by 
standard spectral means. The original treatment by Gel'fand and Yaglom involves 
a limit process from a discretisation approach to functional integration, which, in 
content, is partly equivalent to the remarks here. 

2 There are numerous texts on finite difference equations. An unusual one is Bleich and Melan, 
[f] and a modern one is Elaydi, [2]. The classic work by Atkinson, [3], is a central reference. 
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I treat, at least initially, the simplest set-up that allows me to illustrate the 
essentials. This will be the uniform continuous string of length L vibrating transver- 
sally. An approximation by (equal) mass points takes us back to the precursor of 
Fourier analysis, the subject of countless historical surveys and textbook expla- 
nations. For reference I mention only the classic Rayleigh, [10], and Morse and 
Feshbach, [11]. The modes of this discrete system are, therefore, ancient but I will 
develope them again. Some are given, relevantly, in the basic finite difference text 
by Fort, [12]. There will necessarily be a certain amount of repetition. 

A summary of the discretization, of relevance to the present topic, is given by 
de Verdiere, [13], §9.2. 

In the course of the calculation, I encounter Chebyshev polynomials and deve- 
lope their properties via a two-matrix technique which is convenient, and might be 
novel. 



2. The discrete Gel'fand— Yaglom theorem. 

To make the situation precise, replace the interval [0, L] by z/ + 2 equally spaced 
points, or vertices, two being end, or boundary, points. Label the points by j, 
< j < v + 1 and consider some scalar function, y(j), satisfying either Dirichlet 
(D) or Neumann (N) conditions at the ends, (e.g. [II]), 3 

y(0) = y(u+l) = 0, D 

y(0) = y(l), y(u) = y(u + l), N. 

I discuss the Sturm-Liouville problem which, in its simplest formulation, in- 
volves the eigenvalue recurrence, (e.g. [12]), 

y(j + 1) + (A - V(j) - 2) y(j) + y(j - 1) = , (2) 

subject to boundary conditions, say (1). 

I refer to V(j) as the potential because (2) can be rewritten as the more familar 
looking Laplacian eigenvalue equation, 4 

y(j) = Xy(j). (3) 



For convenience, I will assume that all my functions, eigenfunctions etc. are real, except when 
considering a twisted periodic field later. 
V is the backwards difference operator. 



- 4 VA + V(j) 
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The lattice spacing, h = L/{y + 1), has been introduced by scaling to give a 
'physical' Laplacian and one has the dimensionless quantities, A = h 2 X and V = 

hW. 

The procedure is textbook. Taking D conditions for definiteness, iteration from 
the j = end point, assuming any value of y(l), except zero, yields all the y(j) 
as polynomials in A, in particular, the terminal value, y[y + 1, A). The eigenvalues 
are thus the roots of this polynomial, y{y + 1, A) =0 (e.g. Atkinson, [3]) and the 
determinant (i.e. the product of all the A) of the operator is its constant term, 
y(v + 1,0), up to a factor, which is the essence of the Gel'fand-Yaglom formula 5 . 
The factor involved is unity if the starting term is chosen to be y(l) = 1, as can be 
seen by looking at the A — > oo limit (see later). 

The product of all the physical A is only a scaling factor different and one 
arrives at the discrete Dirichlet result, e.g. [9], 

Vet D = -^y(v + l,0). (4) 

This formula is thus part and parcel of the standard eigenvalue problem. The 
resolvent of (2) is, 

R(\) = ±\ogy(v + l,\), 

with the usual machinery. For example, the sums of the inverse powers of the roots 
follow, a la Euler and Rayleigh, [10] I, p. 279 , as, 

-*(A) = £f + A£-1 + A 2 $:-1 + ... . (5) 

n n n n n n 

The rigorous proof that the discrete formula leads to the original continuous 
one of Gel'fand and Yaglom and of Levit and Smilansky, [6], is given by Forman. 
de Verdiere, [13], also discusses the nature of this limit. 

3. Dirichlet constant potential 

Before continuing to other boundary conditions, I give the simplest application 
of (4) which is when the potential is constant and equivalent to a mass term, fx 2 . I 
then rewrite (2), 

y(j + l)-2cosh2 iy (j) + y(j-l) = 0, (6) 



5 The nature of this constant is where the problem lies in the continuum case. 
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O 9 

where I have set /j, = hp and /j, — A = 4 sinh 7 and which I must solve subject to 
the initial conditions y(0) = 0, y(l) = 1. The roots of the auxiliary equation are, 

m ± = cosh 27 ± sinh 27 = e ±27 , 

and the general solution is, 

y(j) = Am\ + Bm J _ , 

with 

A = -B 1 



2 sinh 27 

This implies that the discrete Gel'fand-Yaglom function (more conventionally called 
the fundamental solution) is, 

sinh27j 

M = linW ' (7) 

evaluated at the terminal point, j = v + 1, which is, perhaps, no surprise in view 
of the textbook continuum analogue. The functions, y(j), written y(j, A), are poly- 
nomials in 4 sinh 2 7 (and hence in A) which can be proved in many ways, one of 
which is the direct iteration of (6). Equation (2) is a recursion formula for these 
polynomials, which are Chebyshev polynomials, as is well known, the definition 
being, 

-r T / _ ^ sinh27(z/ + 1) 

UJcc**,)* >- . (8) 

Pursuing the calculation, the determinant is obtained by setting A = 0, 
^ / n 1 sinh27 n (z/ + 1) n . , 

Det ^> = ^ 5 .nh2 7o ■ " = 2 -^c 0) 

The constant of proportionality is settled by the infinite A limit when the Gel'fand- 
Yaglom function has the explicit behaviour, 

sinh27(. + l) ^ ogh 

sinh 27 v u\) 

The eigenvalues themselves are determined by, 

sinh27(z/ + 1) = 

or 

mri 

7 = 7n = 257+I) ' 
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and so 

— _ 2 4 2 nn 

K = » + p sin 2 + , n = l,...,z/, (10) 

which is the textbook result, e.g. [12]. Equating the determinant (9) to f| n gi ves 
a standard product formula, e.g. Bromwich, [14], p. 211, which comes up later in §13. 
Furthermore, the sums of inverse powers of the roots, (5), yields finite summations 
for powers of cosecants, e.g., typically, 

X>- 2 ^^ 2 -i), ( n ) 

n=l 

which are very old and are simple examples of a wide class of trigonometric summa- 
tions obtainable in many ways. 6 As p — > oo (the continuum limit) this sum becomes 
Euler's result, ( R (2) = n 2 /6. 

The eigenfunctions follow by noting that the fundamental solution, y(j, A), (7), 
satisfies the equation (6) with A = A n and obeys the Dirichlet conditions. The 
eigenfunctions are therefore, 

2/nU) = sm ^py 5 n = l,...,z/, 

of which there are v, this being the number of 'dynamical' points. We therefore 
reach the standard mode properties, e.g. Fort, [12], Spiegel, [17]. This route is not 
a novel one. 



4. Neumann conditions 



As a warm-up for the Robin case, I consider Neumann boundary conditions 
(1) which can be written, Ay(0) = 0, Ay(v) = 0. If V(j) were uniform, (2) would 
be satisfied by Ay(j) and the problem translated into a Dirichlet one (see §8) but, 
because of the j dependence, this is not possible and it is necessary to treat (2) and 
its difference together. 7 This is most neatly expressed using 2x2 matrices as in 
[3], [9] and elsewhere, e.g. [7], although a little differently. It might be considered a 
'phase space' representation. 

Defining, 

^ = L- +D ) ■ <i2) 

6 I attempted a few comments and gave some references in [15]. See also Berndt and Yeap, [16]. 

7 This is a common device in the theory of ordinary differential equations. For difference equations 
see e.g. Porter, [18], Goldberg, [19], p.233 Ex.4., Elaydi, [2]. 
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the recurrence under study is the first order one, 

T(j)-M(j)T(j-1)=0, (13) 

where 



M(j) = 



(14) 



1 

-1 V(j) + 2-X 
1 
-1 2cosh27 j 

whereby 7 • is defined. One of the equations is just an identity. Note that det M(j) = 
1. 

The Sturm-Liouville Neumann boundary condition is Ay(0) = 0, Ay(v) = 
which means, choosing a normalisation, 

T(0)=Q, T(i/)<x(j). (15) 

The eigenvalue procedure is to iterate (13) up to T(z/) starting from, T(0) i.e. 

T(i/) = M(i/)M(i/-1)...M(1)T(0), (16) 

and then impose the condition, (15), on T(z/). (Remember, the Ms are functions 
of A.) 

The roots of the polynomial, 

P(u,X) = (l, -1 ) Y(i/) = Ay(v, A) , 

are then the eigenvalues and the determinant is P(v, 0), which is the required the- 
orem, [9], 

De W = ^;A 2 /( i /,0). 



5. Robin boundary conditions 

Having treated pure Neumann, it is not much more difficult to sort out Robin 
conditions, called General Local in [9]. The recurrence is still (13) but with the 
boundary conditions, 

T(0)=( 1 |J=T inJ TH^ + ^eT^, (17) 
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where a and (3 are the Robin parameters defined by, 

Ay(0) = a y{0) , Ay(u) = -(3 y(v + 1) . 

The first condition in (17) is chosen, and the second is imposed. 
Defining an 'adjoint', 

Tt(j) = T0V, 

in terms of the symplectic metric, J = ^ ^ , the eigenvalue polynomial is, 
therefore, the matrix element, 

(-1, i + p)y in (v)=rt ut y in (is) 

=TL%)I(,-l)...M(l)T in) 

and the product of the A eigenvalues is proportional to T\ ut T in (v) , evaluated at 
A = 0, the constant of proportionality being the inverse of the coefficient of the 
highest power of A in (18). For very large A, M, (14), approximates to, 

M{j)~(° n \) , (19) 



.0 -A, 

when the right-hand side of (18) becomes (1 + /3)(1 + a) (—A) 1 ' so the determinant 
is, after scaling to the physical eigenvalues, A, 

Det*(5,£) = ^7 (1 + /3) 1 (1 + Q) TL f[M (j)T in , (20) 

where M is the matrix M evaluated at A = 0. For future use I have also introduced 
the physical constants, a = a/h, (3 = /3/h. 

I note that the symplectic product is just the Casoratian (discrete Wronskian), 

[2], ' , ' 

tI(j)t 2 0") = t 1 ovt 2 0"), (2i) 

of two solutions, T-^ and T 2 , of (13), which is a symplectic development because, 

MJM = J, 

and so (21) is uniform, i.e. independent of j. This is a neater proof of this fact 
than the usual ones, e.g. Fort, [12]. 

The equivalence of a 2 x 2 real matrix formulation and a three-term recurrence 
relation is expounded by Atkinson, [3], §3.5 involving a geometrical interpretation 
of symplectic action. 

I take up the general formalism again in §10 and turn now to an elementary 

case. 
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6. Constant potential 



As the simplest example, I again take that of a constant potential. Then M(j) 

is, 

M(j) = M = ( 1 J , (22) 

yJ> V-l 2cosh2 7 ; ' 1 ' 

with 7 as before, i.e. 4sinh 2 7 = [i 2 — A. 

It is shown in Appendix 2 that the power M v is given by, 



M v = 



-U v -! U„ 



in terms of Chebyshev polynomials, (8), of argument cosh 27. Then, 



TLM"T 4n = (-l,l + /3)( Vv ~ 2 



u v -A f 1 

U v J \l + aj (23) 
= (a + + ap) U v - {a + - X)U l/ _ 1 , 



which is to be substituted into (20), after setting A = to give the determinant. The 
formula is symmetric under interchange of a and /3, as it should be by geometric 
symmetry. 

As another check, the Dirichlet choice, 

reproduces (9). 

Incidentally, if M were more general, rather than diagonalisation, it would be 
easier, for iteration purposes, to set, 

M" = al + bM, (25) 

and compute a and b in terms of the eigenvalues of M, e.g. Goldberg, [19]. 
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7. The characteristic polynomial and Euler-Rayleigh sums 



The Euler-Rayleigh sums, analogous to (11), arise from the expansion of (23) 
in powers of A, i.e. of — 4sinh 2 7 (for simplicity I set [i to zero) which is easily 
accomplished by, say, using Bromwich, [14] chap. IX or the relation to Chebychev 
polynomials. I find, 



1 out 1V1 1 in 



= w+«+«£^( 2 7>Ar-. (26) 

+ (a? - i(a/3 + a + + 2)a) £ _ * " 1 ) 1 . 

The v eigenvalues, A n , are the roots of this characteristic polynomial and it is next 
required to expand its logarithm, which, for low powers, can be done directly by 
hand. As the simplest case I give 

V cosec^J = 2 3 " 2{aP + " + g + - ^ + M f + " + g + 2) , (27) 

3((l + ^)a/3 + a + /3) 

where I have set 

K = 4 sin2 ^n, 

which defines 6> n . 

The continuum limit h — > 0, z/ — > oo is not without interest and is discussed in 
a section 9. 



8. Neumann conditions revisited 



As the free N case is not given in Fort, [12], I give, for pedagogic completeness, 
the conventional calculation by noting, first, that the N conditions, (1), can be 
written Ay(0) = Ay(is) = 0. So, defining = Ay(j), one has, from (3), 

-ivA^) = A^) (28) 

with (p(0) = (piis) = 0, which is a D problem on v + 1 vertices but with the original 
spacing, h. The D-eigenf unctions are 

00) = sin— — , n=l,...,i/-l, 
9 



and hence the iV-eigenf unctions are,* 

y(j) = A" 1 cj){j) = A~ L sin 



nn(2j - 1) (29) 

~ cos , n = 0, 1, . . . , v — 1 , 

2v 

up to a numerical factor and possible additional constant. The eigenvalues are, 

A n = -^sin 2 ^, n = 0,...,z/-l. (30) 

Again, there are v modes, including the uniform zero one, n = 0, which cor- 
resonds to a constant of integration in (29). [n = v gives a vanishing mode.) 

Before going on, it would be best to check, if the pure Neumann 

determinant, for the free case with mass, agrees with the above mode structure and 
the Robin expression, (23). Effectively I start again. The initial condition that fixes 
the Gel'fand-Yaglom function is z(0) = 1 and z(l) = 1, i.e. Az(0) = 0. (This is 
Forman's z.) The general solution is again, 

z (j) = Ae 2lJ + Be~ 2 ^ , 

and the conditions imply, 

A + B = 1 

Ae 27 + Se" 27 = 1 , 



which solve to, 



so that, 



A = , B = 



2 cosh 7 ' 2 cosh 7 ' 



z 0) = = V j-i = W j-i 

Az(j) = 4 sinh 2 7 "^ j = AVUj_ 1 . 



(31) 



cosh(2j — 1)7 
cosh 7 
2 _ sinh 27 j 
sinh 27 

where Vj is a Chebyshev polynomial of the third kind, [21], and all Chebyshev argu- 
ments are cosh 27. Appendix 1 contains some relations for Chebyshev polynomials 
couched in the two-matrix language and (31) can be obtained more rapidly using 
this. 



8 If you use Jordan, [20], be aware that there is an error on p. 117 that is carried forward. For 
example, on p. 124 the sum of cos(ie + b)(j> is incorrect. The upper limit should be n — 1. 



10 



Applying the eigenvalue restriction, Az(v, A) = Az{y) = 0, yields the condition 
7 = 7 n = nrri/v, n = 0, 1, . . . v — 1, and the eigenvalues are, 

4 7TTZ 

K = + -^2 sin2 ^7 ' n = °> • • • ' v ~ 1 ' ( 32 ) 

consistent with (30). The eigenf unctions, (29), also follow trivially from (31). 

One sees from (31) that the eigenvalue condition is the same as the Dirichlet 
one, except for the replacement v — > v — 1 and for the factor 4 sinh 7 = /z — A. 
This factor is responsible for the n = mode which, in the massless case, is a zero 
mode. 

The Neumann determinant is then, 9 

De W(/«) = pjM^o), 

where the numerical factor follows on the limit, 

4 sinh 2 7^^ ^4 sinh 2 7(2 cosh 2 7 ) ~ (-A)" . 
sinh 27 

Equating the two forms of the determinant yields the same product formula as in 
the D case. 

The determinant also agrees with the Robin formula, from (20), for a = (3 = 0. 
(This is, of course, simply a check of algebra.) 

9. The continuum limit 

Comparisons with known results can also be obtained by considering the con- 
tinuum limit, an historical motivation for discretisation. Again as an example, I 
consider the Robin determinant (20) with (23) in the limit h — > 0. To get the 
leading divergence, the lowest power of h is required in the expression multiplying 
l/h 2u . As h — > one has the limiting behaviours, 

2 sinh 7 ~ 27 ~ /z = hjl , 2^ u ~ tifiv ~ JiL , 

and therefore by inspection of (23), one sees that the leading term is of order h. 
Extracting this gives, 

T^" M v T 

— — _ 2 ,sinh7lL 

(33) 

— > (a + p ) cosh [iL + [ap + /i ) — — — , 



This appears to differ by a factor of 1/4 from Forman's formula, [9]. 



11 



which agrees with an expression in [22] for the continuum case. 

Related is the limit of the simplest Euler-Rayleigh eigenvalue sum, (27). Re- 
verting to physical quantities, 

h 2 V — = V — -> 3 (^ + ^) + ^ + 6 / 34 n 
which is also given in [22]. 



10. Non— uniform potential. The propagator 

Difference equation Sturm-Liouville theory is well developed and can be pur- 
sued by analogy to the continuum version, e.g. Fort, [12], Levy and Baggott, [23]. 
In fact Sturm obtained many continuum results via a discrete route, although this 
was never published. 

In this section I wish to develope and summarize the previous matrix formula- 
tion, see (12), (13), (22). I consider (13) as a Schrodinger equation for a two-state 
system with a discrete time labelled by j, and rewrite it by defining a matrix 'prop- 
agator' K(X;j,f), 

j 

K(X;j,j') = e(j,j') M(k) (35) 

k=j'+i 

as 

r(j) = K(\;jj')r(j>), j> 3 ' 7 (36) 

which propagates forwards from j' to j and acts as a transfer 2x2 matrix. In the 
simplest case, the matrix M is given by (22). The form, (35), is an equivalent of 
the time-ordered exponential solution in time-dependent perturbation theory, but 
here 'vertex-ordered'. The propagator, K(X;j, 0) is sometimes referred to as the 
state transition matrix. I denote it by K(X;j). The basic theory is given by Elaydi, 
[2] §3.2, but my treatment is modified a little and also deals, particularly, with a 
symplectic invariant propagation. 

For consistency, the initial condition, (i.e. the first empty product in (35)), 

K(X;j,j) = l, (37) 

has to be taken. The step function ensures that K(j,j') = for j < j', corre- 
sponding to causal propagation. The semi-group property, 

K(X;j,f)K(X;f,f) = K(X;j,f), 
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(no sum on j') and symplectic invariance, 

K(X;j,j')JK(X;j,j') = J, (38) 

also hold. 

K satisfies the equation of motion, 

K(X;j,j') = EK(X;j - 1,/) = I6 jtj , + M(j)K(X;j - 1,/) , (39) 

where the first term arises from the 9 factor in (35). A matrix which satisfies (39) 
is a fundamental matrix. 

Iteration of (39) gives a power series expansion, 

K(X-jj') = lS jd , + M(j)5 JtJ , +1 + M(j)M(j - l)S jJ/+2 + . . . , (40) 

which is quite equivalent to (35). It also follows from the decomposition, 

0(3 J') = S jd , + 6 jJI+1 + 6 jtjl+2 + ... 

= A 6 J,J" 

obvious graphically, arithmetically and in (v + 1) x (y + 1) matrix form. It is the 
discrete version of the distributional operator statement that the ^-function is the 
integral of the 5-function. 

If M(j) is constant, then, trivially K(X;j,j') = 6(jJ')M j - j \ either from (35) 
or read off from (40). 

If an 'unperturbed' propagator, K Q , is defined by, 

KqUJ') = 16 jtj , + M (j)K (j - 1,/) , (42) 

then, 

K(X;j,j") = K (j,f) + K (j,f)(M(f) - M (f))K(X; f - lj") , 

where j' is summed over from 1 to z/, can be considered as a perturbation expansion. 
If M Q is constant, 

K(X;j,j") = Mt j " + Mt*' (M(j') - M )K(X;j' - 1,/') . (43) 

The propagator, K(X;j,j') is defined independently of any boundary condi- 
tions which are incorporated, in my approach, by constructing the symplectic scalar 
products, 

PW = ^L(j)K(X;^j')T in (f) 
13 



These are polynomials in A and, because of the uniformity of the Casoratian, 
are independent of j. The boundary conditions are given by T in (0) = T in and 
T out( iy ) = T out: as given in (17). T in (j) is the solution of (36) for the 'in' condition 
and T out (j) that for the 'out' one. 

The vanishing of P(X) determines the v eigenvalues, A n , (n = 0, — 1). 

This characteristic polynomial reads, in the extreme cases, 

p(a) = tLm ^(i/) = tL(o) T in (0). 

The A dependence is contained in T in (v) or in T| ut (0). 



The full determinant is the normalised -P(O), 
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Det = t P '°' . (45) 
t' at 

1 out^ 1 in 

All this we have had before in particular cases. 

To expose the parameter A, and to enlarge on the formalism, it is helpful to 
split the driving matrix M as, 

M(j) = B(j)-\A(j), (46) 

where, 

/ 1 \ /o o\ 

B(j) = , A(j) = A =[ , (47) 

\-l V(j) + 2j V° V 

with 

BJB = J, AJA = 0, AJB = -A. (48) 

Then consider two fundamental matrices, K(X;j) and K(fx;j), and make the 
usual construction, cf [3], 

K(»;j + 1) JK(X; j + 1) - K(fi;j) JK(X;j) 

= K(»;j)((B(j) - nA) J (B(j) - XA) - J)K(X;j) 
= (X-fi)K(^j)AK(X;j). 

Summing over j from to v — 1 (i.e. performing the inverse A -1 ), one gets 

v-l 

K(p; v) JK(X; u) - J = (A - //) ^ K(p;j) AK(X;j) . (49) 

3=0 
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P(X) is the analogue of an S'-matrix element. 
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In this equation, A and \i are any two parameters. I now restrict them to being 
eigenvalues, that is, solutions of the polynomial equation -P(A) = 0, or, 

T L0') T m0') = 0> V J, 

which implies that the 'out' eigenvector T out (j) is the same (up to a constant 
factor) as the 'in' one, T in (j), for each eigenvalue 11 and I can denote both of them 
by T (j, A). In particular, T out (0) and T in (u) are independent of the eigenvalue. 

In this case, mutiplying (49) by the boundary (eigenvalue independent) vectors 
T in , on the right, and T out on the left, the left-hand side vanishes, 

= T out (0, fi) J r in (v, A) - f out (u) J T in (0) 
= T(0)JT(i/)-T(i/)JT(0) 
= 

and so, 

u-l 

(A-rf^T(^)AT(j,A) = 0, 

3=0 

with the usual conclusion that eigenvectors with different eigenvalues are orthogonal, 

^T(j, / i)iT( J ,A) = 0, ^A. 

3=0 

In the circumstances of the present paper, the matrix A is the projection onto 
the lower components of T(j), i.e. onto y(j + 1), and so I have regained the usual 
orthogonality, 

V 

Yl 2/0"> K) VU> X m) = Pn S nm » 0<7l,TO<Z/-l, (50) 

where p n is a normalisation. 

A standard procedure then yields the completeness (or dual orthogonality) 
relation, 

y& K) y(f, K)pn x = $jj> , 1 < h i' < v . 

n=0 



11 Equivalently, the eigenvalues are simple. 
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The appearance in (50) of a sum over the vertices, j, leads us to the tradi- 
tional matrix approach, (e.g. Rayleigh, [10], Atkinson, [3], Chap. 6, Gantmacher 
and Krein, [24] ), which takes the entire set of (dynamic) values, y(j) (1 < j < v) 
as the components of a column z/-vector, y, and writes the collection of difference 
equations, (2), as a v x v matrix equation of the familiar eigenproblem form, 

L y = Xy . 

The polynomials, y(j, A), are then related to the Jacobi determinant det (L — Al). 
I extend my present formalism a little to reflect this perspective, which has already 
shown up in (40) and (41). 

It is sometimes convenient to employ the operator formalism as in finite di- 
mensional quantum mechanics, due to Schwinger and Weyl, [25], cf Floratos [26], 
and set e.g., 

(j\K\j') = K(j,j'). 
I retain K as a 2 x 2 matrix in phase space. Then the recurrence is, 

T = MX, 

with M a subdiagonal matrix, 

(j\M\j') = M(j)5 jJI+1 , 
and the series (40) translates into the simple operator equation, 

K = 1 + MK, 

or, formally, 

1 - M 1-B-XA 
The elements (which are matrices) of the powers of M are correctly vertex ordered. 

I will not pursue this formulation any further at this time except to say that the 
stepping matrix has just ones along the subdiagonal and represents the translation 
operator, often denoted by E in finite difference calculus. The Heaviside matrix, 
©, having 9(j,j') as elements, is triangular with ones in the left-hand part, and on 
the diagonal. It is related to E by EQ = — 1. 
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11. The 6 potential on the interval 



A very basic example of a variable potential is one that is non-zero at only one 
vertex, i.e. V(j) = vSj k ,. Then, in the product form, (35), of the propagator, only 
one term k = k' will be different from the rest. The remaining products (powers) 
can be dealt with as before, in §6, and an explicit expression found for the transition 
operator K(X; v), say. 

As this is just meant for illustrative purposes, I choose a value of k', namely 
k' = 2, that results in a simple formula. 

In this case, for Dirichlet conditions with (24), I find the polynomial in A, 

xL^(A; v)T in = U v (l - A/2) + v{2 - X)U V _ 2 (1 - A/2) 

= U v (l - A/2) + v(U v _ 1 (l - A/2) + tf„_ 3 (l - A/2)) , 

in terms of Chebychev polynomials (the 'unperturbed' functions), (8) with 2 cosh 27 = 
2 — A. The eigenvalues are easily determined numerically and, for a small number 
of vertex points, even analytically as functions of the strength of the potential. 

Equation (51) is proved and explored in Appendix 2 where it is extended to 
the case when all the Vj are populated. 

As particular quantities, the sums of the inverse eigenvalue powers can again be 
computed by expanding the logarithm of this polynomial, which, to lowest orders 
is, 



ri ut K{\-v)T in = (v + l)+2v(v-l)-^(v(v+2)+2v(v 2 -3v+3)y+. 
on using, 

Encash 2 7 ) = (u + 1) ( 1 - \v{y + 2) A + ... J . 



(52) 



One then finds, exact in v, 

^ 1 _ u + 1 (is(is + 2)+2v(is 2 -3is + 3)) 
^ n K~ ~ 6 v + l + 2v(v-l) ' 

n=0 n v ' 

which generalises (11). 

This identity is an example of a general class of identities discussed in the 
interesting work by Annaby and Asharabi, [27], where other references can be found. 

The determinant is just the constant term in the polynomial, (52), 

Det D = 1/ + 1 + 2u(i/ - 1) , 
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and a zero mode, A = 0, occurs when v = —{y + l)/2(z/ — 1). When v takes the 
same value with the opposite sign, the final eigenvalue, A„_ 1; equals 4. 

I remark that perturbation theory on the lattice has been considered by Actor 
et a/, [4]. 



12. The vacuum energy 

A rather different technical eigenvalue problem is the calculation of the Casimir 
energy and I present a quick treatment as a simple and explicit use of the eigenvalues. 

The Dirichlet vacuum energy of a free scalar field on T x I v can be evaluated 
in closed form as, 



2 

A 

1 \ . nn 



sin 



^Sr" 2 ^ + 1 ) (53) 

= - ( cot — — — - - I 
h\ 4(i/ + 1) 

21 1 71 

If one views the lattice calculation as a regularisation of the continuum one, 
the — n/2AL term is recognised as the value given by the (^-function technique while 
the first two, ultimately divergent terms, being non-universal, dependent on the 
regularisation, should be discarded in some way, if one is concerned just with the 
interval, [0, L] on its own. 

The paper [4] contains a full discussion of the expression (53) and I will not enter 
into any more details. This reference also contains other arrangements, including a 
discrete version of the Casimir piston. 

The Neumann energy is, likewise, 

1 ^ . wn 

n=l 

(54) 

h v Au ' 

2L 2 1 TT - „ 

+ . . . , h->0. 



nh 2 nh 2h 2AL 
The other boundary condition usually considered is the periodic one. This is 
given in Fort, [12], Chap. XV. It is convenient, this time, to arrange the v points, 
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< j < v — 1 on the unit circle and impose the periodicity conditions y{y) = y(O), 
y(— I) = y{v — 1) which relate values outside the proper range of j to those inside. 

The analysis is slightly different depending on whether v is even, v = 2k + 2, 
or odd, v = 2k + 1. In both cases there are degenerate modes, cos(2nirj jv) and 
sin(2n7rj/z/), for < n < k with eigenvalues, 

— 4 2 7m 
A = — 5- sin — , 
br v 

where the gap, h = 2%/u. The value n = gives the one uniform zero mode. If u 
is even, the single mode cosnj must also be added. (This alternates between plus 
and minus one as the points around the circle are traversed and corresponds to a 
wave of infinite frequency in the continuum limit.) The total number of modes is 
always v. 

In exactly the same way as above, the vacuum energies are, 

2& + 2/A nn 1 
E ^ = —{^ Sm 2kT2- + 2 

x n=l 

2k + 2 tv 
cot 



2tt 2(2/c + 2) 

2k + 1 nn 

-C/oi.,1 = > sin — 

2fc+i n 2k +1 

2k + 1 7T 

cot 



or 



2n 2(2Jfe+l) 

S p = -oot--> )? -- + ... > (55) 

in both cases, as expected. Again, one sees the continuum zeta value of ( R (— 1/2) = 
— 1/12 appearing as h tends to zero. 12 

Fort, [12], also discusses anti-periodic (real) functions. However I will be a 
little more general and analyse a system that, in the continuous limit, amounts to 
an Aharonov-Bohm flux running through the circle. This is mimicked by imposing a 
phase change on circulating the flux and leaving the equations of motion unchanged. 



12 There is a puzzle here. In the continuous case the periodic modes on a circle are the union 
of Dirichlet and Neumann modes on an interval of size half the circumference. One might, 
therefore, expect to see evidence of this, even in the discrete case, as h — > 0. In fact this works 
for the terms of order h~ 2 and h° in (53), (54) and (55) but not for those of order h —1 . In 
order for it to work, the relevant term in (54) should read just l/2h to cancel that in (53), on 
addition, to give (55) but I could not achieve this. 
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In the quantum case, the wave function is complex and exponential functions 
are very convenient. 13 I therefore consider a function, defined on the points, 

j, and satisfying the twisted periodicity condition, 



if,(v) = e 2 ">(0) , tp{u- 1) = e 27Tia i/j(-l) . (56) 
The modes on the discrete circle are, 

^(j) = e ^i(n+a)j/u^ re = o,...,i/-l, 0<a<l, 

with corresponding eigenvalues (h = 2tt/v), 

- 4. 2 7r(n + a) 

A = — sin 2 -i '- , 57 

h v 

and vacuum energy (with a factor of two from the complexification) , 

. « 2 ^-^ n(n + a) 

E(a) = - > sin -i '- 

n=0 

2 h h (58) 

= — cosec — cos — (2a — 1 
h 4 4 

8 1 

E(a) must be extended beyond a = 1 using periodicity, i.e. i?(l + a) = 

The constant term agrees with the result (the periodic Bernoulli polynomial, 
B 2 ) that arises in the continuous circle limit, [28]. When a = one regains twice 
the real periodic value (55). It might be of interest to give the full formal expansion, 



m=0 

which I have not seen elsewhere. 



^> = 2 E^ 2m(a) (!) 



2m-2 



13 It is, of course, possible to retain a real description by doubling up the fibre to an SO (2) one. 
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13. Direct determination of determinants 



It is helpful to have specific values for comparison or limit purposes and I 
proceed to evaluate the determinants of the free systems directly from the eigenval- 
ues which have just been used. I also look at the continuum limit and some zeta 
regularised values. 

An organisational point is perhaps required. In deriving the answers, I have 

used known expressions for some infinite products (e.g. Bromwich, [14]), which could 

be obtained, according to the results of the previous sections, from the determinants 

and so these evaluations might be considered superfluous, at one level. However I 

present them as independent checks. I also include the twisted periodic values. 

To repeat, Dalembert's equation is, 

1 



~VAy(j)+JI 2 y(j) -Xy(j) = 0. 



(59) 



For the interval, I will again use the real form of the eigenfunctions and the 
calculation is immediate. 

The D-determinant on the L-interval using the eigenvalues (10) is, 



Bet D (fj,) 



2v v 



n 



1 



sin 



2(1/ +1) 



2v r 2v+V 



sin 



n 

n=l 

1 sinh(z/ + 1)27 



irn 



2(1/ + 1) + 4^ 



fj, = hji 
1/2 



(60) 



h 2u sinh 27 
The dash on the second product means that the n 
and I have set \x = 2 sinh 7. 
The massless values are, 

1 



v + 1 term is to be excluded 



Det D (0) 



h 2 » 
1 



li 2 »- 1 2 2L ' 



The N-determinant is, using (30). 
Detjv(7Z) = 



n 



irn 



n=0 
2v v 2v-\' 



sm 1 — a' 



n = l 



n ( . 2 nn 1 
Sm 217 + 4 



/// 



1/2 



(61) 



-7— 2 tanh 7 sinh 27Z/ . 

h 2u 
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In the massless limit the determinant vanishes and it is conventional to remove the 
offending zero mode giving the modified determinant, 

Before discussing these results, I give the twisted periodic expressions. 
From (57) (h = 2tt/v), 

V / n =0 ^ 

= ^7(cosh2 7 z/-cos27ra) (62) 
Detp /2 (a,0) = -^4 sin 2 yra , 0<a<l, 
Det / P (0,0) = ^^4(27r) 2 , a = 0. 

For the a = case, I have removed the complexification squaring and, for extra 
generality, I have included the mass term. 

For comparison, some determinants, computed from the bare (^-function regu- 
larisation, are well known to be, 

^ . . sinhuL 

Det' CJV (0) = 2L 
DetJ P 2 (ct,0) = 4 sin 2 yra 
Det' CP (0,0) = (2tt) 2 , 

and one sees that the lattice determinants are proportional to the zeta values, in 
the continuous limit. In particular, 

DeUp) = 1 ^ + 1)2, 
mPV h 2v sinh27 

1 sinhuL 

-» To-. n — » h^O 

h 2u+1 -p 

Forman introduces a twisted periodic condition, denoted B 5 in [9], which has 
a complex multiplying factor 6 instead of e 27Tlct in (56) and Theorem 2.6 gives its 
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determinant. Evaluating [9] equation, (2.23), in the free case I find, 



1 6(1 -8) : 

Det, - 



J B 5 



h 2 » (1 + \5\ 2 ) 



which vanishes when 5 = 1, as a check, but does not agree, apart from the lattice 
scaling factor, with (62) when 6 = e 27riQ . 14 

14. Conclusion 

A good deal of this paper is expository but there are some novelties. It has 
been emphasised that the Gel'fand-Yaglom formula for the determinant in the dis- 
crete case is a standard component of Sturm-Liouville and orthogonal polynomial 
theory. I have rewritten this in a neat 2x2 symplectic matrix formulation, slightly 
different from the usual one, and have calculated the determinant for Robin bound- 
ary conditions for a constant potential, as a basic example. The continuum limits 
have been discussed in a simple-minded way and comparisons made with the work 
of Forman, [9], revealing some minor discrepancies. For Dirichlet conditions, the 
determinant for a 5 potential was evaluated exactly, highlighting the significance of 
Chebychev polynomials as 'unperturbed' Sturm-Liouville solutions which is further 
explored in the Appendices and which contain some technical advances. 

The calculations could be broadened to include the general Sturm-Liouville 
operator, and higher order equations. 

Added note 

The Green functions for free propagation on the discrete interval (or 'path'), for 
various boundary conditions, have been obtained in terms of Chebyshev polynomials 
by Bass, [29], Chung and Yau, [30] and Bendito, Encinas and Carmona, [31]. 



14 Curiously, if one of the 5s is replaced by 5 1 , then agreement is found, apart from a factor of 
two, which seems too much of a coincidence. 
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Appendix 1. Chebyshev Polynomials 



I give some basic results for Chebyshev polynomials in a way that reflects the 
procedures of this paper. These polynomials occur in the Chebyshev-Gauss scheme 
for mechanical quadratures. Probably the most economical way of defining them is 
through the recursion, 

P n+1 (x) - 2xP n (x) + P^x) = , (63) 



or 



VAP n (x) = -(2-2x)P n (x), (64) 



the different kinds being selected by the 'initial' values, P and P x , e.g. [21]. 

As a slight novelty, I use the matrix description adopted in the main body of 
this paper. So, introducing the two-vector, 



TT f X _ ( P n (x) 



the three-term recursion (63) becomes the two-term matrix one, 

U n (x) = C(x)U n _ 1 (x), 

i.e. 

U n (x) = C n (x)U (x), 

where 15 

°w=(! D- (6B) 

This allows me to specify the initial conditions neatly. For example, the 'Dirichlet' 
vector, 

n D (^(°' 

produces Chebyshev polynomials of the second kind, P n (x) = U n (x) while the 
Neumann one, 

gives third kind polynomials, P n (x) = V n (x), [21], [32]. 



15 C is what I call M in the main body of this paper. 
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Many relations between the various kinds can be obtained by combining initial 
conditions. 

The first iteration of the initial vector ( I gives ( I i.e. minus the 

1 

Dirichlet one and therefore yields —U n _ 1 after n iterations. So, writing I 

+ ^q^> one obtains the relation V n = U n — U n _ 1 = AU n _ 1 . From which, for 
example, AV n = VAU n , helpfully relating Neumann and Dirichlet. 
The Robin choice, 

likewise generates the D-N combination V n + aU n = (1 + a)U n — U n _ 1 . 

The different polynomial combinations are encapsulated in the form of the 
power, C n , 

(U_ l = 0) which can be written, 

C n = U n A + U n _ 1 J+U n _ 2 A', (67) 

with A' = A — 1. From (66), by taking the determinant, one finds the Chebyshev 
identity, 

U 2 n _ 1 -U n U n _ 2 = l>Q, 

which is a (known) statement of a Turan inequality. 

From the Z group composition rule, C m C n = C rn+n 1 the combination relation, 

Um+n = U m U n — U m _ 1 U n _ 1 , 

can be deduced and the SU(2) character Clebsch-Gordan series, 



u m u n = J2 u k> 



k=\m—n\ 

readily follows therefrom by iteration. All these relations have trigonometric deriva- 
tions. 
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The generating function can also be transcribed as follows. The matrix equa- 
tion, 

'l -t 



(l-ct)- 1 



t l-2tx 

1 / 1 - 2tx t 



l-2tx + t 2 \ -t I 

oo 
n=0 

and (66), yield the normal generating function for U n . (Incidentally, to check the 
top left entry, one has to use U_ 2 = — 1, U_ x = 0.) Other identities can be deduced 
in a similar fashion. 



Appendix 2. Calculation of matrix elements 

I give some details of the computation of the matrix element polynomial, 
-P(A) = Y\ ut K(\;v)Y in , whose vanishing determines the eigenvalues. 
I write, in dyadic form, 

P(A) = Tr(KT in ®TL)=Tr(KQ), (68) 

where the matrix Q takes the specific forms in the D and N cases, 

Qd= (° j)=a, 0w= (-j ;). 

This is a trace formula for the Gel'fand-Yaglom function. Other expressions 
are determinant ones, [3], p. 95, and also, e.g., [9], [7]. 

The general form of K(X; v) is a product of v matrices, (35), but as a first 
example I treat the case discussed in §11 where all matrices are identical, except 
one. I use the notation of the previous Appendix and write K as, 

K{\-v) =C m (x)C(y)C n (x) , m + n + l = v, 

where C(x) and C(y) are given by (65) with 2x = 2 — A and 2y = v + 2 — A. Then, 

K = C m (x) C(y) C n (x) = C v {x) + C m (x) (C(y) - C(x)) C n (x) 
= C v {x) + vC m {x)AC n (x) 
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and hence for (68) in the case of D conditions, 



P D (X) = Tr(C"(x) A) + vTt (C m (x)AC n (x 
= U iy + vU n U m 
= U u + v{U n+m + U n _ 1 U m _ 1 ) 
= U v + v(U v _ 1 + U n _ 1 U v _ n _ 2 ), 



)4 



(69) 



where I have used (66). Setting n = 1 gives the result (51). 

Equation (69) gives the first terms of a perturbation expansion, to elucidate 
the general nature of which, in a direct way, I consider the case when there are 
two distinguished matrices in the product (35), and a more systematic notation is 
required. 

The vertices of the interval have been labelled by j which runs from to v + 1 
with the vertices 1 to v being dynamic. 16 In this range, let the vertices j 1 and j 2 
be singled out to correspond to matrices C{y^) and C{y^). Then, taking j 1 > j 2 , 
and noting that C(yj ) = C + Vj A, with C = C(x), I find, 



I have added the post factor of A to bring out the fact that it is the combination 
C n A that enters, and this takes the form, 



The final act of taking the trace picks out just the lower right corner term and 



K(\- v)A = C-ii C{ Vji ) C^-'CiyJ C^A 

= C v ~^ (C + v H A) C^-^-^C + v j2 ) C 3 *- l A 
= C 1/ A + v 1 C u ~^ AC^' 1 A + v, C-i'AC'^A 
+ u, v; C V ~^ AC^-i*- 1 AC 3 *' 1 A. 

Jl J2 




which is preserved under multiplication, 




so, 



Tr (K(X; u)A) = U„ + U V _ H Vji U jx _ x + U V _ H v ja U^_ x 

+ U »-h V h U h-i*-l V 3* U h-l' 



(70) 



In graph theory language these would be internal vertices. 
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The structure when more vertices are marked is clear and the general expression 

is, 

V V 

Tr(K(\ ] v)A) = U v +Y,U v - jl v ji U ji _ 1 + £ U^^U^^U^ 

ii =1 ii>i2 =1 

ii>i2>i3 =1 
+ ... 

(71) 

This equation also follows (equivalently) from the complete iteration of the 
conventional-looking (43), which I rewrite here, in the notation of these Appendices, 

V 

K(X;j) = &{x) + C j - j i(x)A Vji K{\- 3l - 1, 0) , 2x = 2 - A. (72) 
Ji=i 

This can be given the usual propagation interpretation of a perturbation series (al- 
though finite and exact) with the Chebyshev polynomials acting as free propagators. 
An obvious graphical representation can be set up. 
Setting A to zero in (71) gives the determinant, 



J2 



Det D = v + 1 + ( u ~ k + 1) k v dl + ( u ~ k + !) (k ~ k) k v h 

j'i=i j'i=ij 2 = 1 

+ Y1 Yl ( u ~ k + !) Ch - J2) O2 - k) k v j 1 v j 2 V J3 

ii =1 i2 =1 i3 =1 

+ ... 

+ v 1 v 2 ...v v . 

(73) 

The upper limits have been extended to the diagonal values using the vanishing of 
the summands there, but I make no use of this at the present time. 

As a numerical illustration, the characteristic polynomial when v = 3 is, 

P(\) = -A 3 + A 2 (6 + S t ) - A(10 + 4S t + S 2 ) + 4 + v 2 + 3S 1 + 2S 2 + S 3 , 

where 

S 1 = v 1 +v 2 + v 3 , S 2 = v x v 2 + v x v 3 + v 2 v 3 , S 3 =v 1 v 2 v 3 . 
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It is interesting to note that, for a symmetric potential (v 1 = v 3 ), P(X) has the 
linear factor, (A — — 2), which is related, presumably, to Borg's reconstruction 
theorem. 

Neumann boundary conditions can be handled in a like manner. The required 
polynomial turns out to be 



Tr (K(X; u) Q N ) = AV,_, + £ V v _ h v jx V jx _ x 

Ji=i 

V 

+ y v„_ i Ua ,n, Va _i 

v Jl Jl Jl J2 1 J2 J2 1 ('7/| N l 

ji>j 2 =1 ' ^ 

V 

j'i>J2>J3 =1 
+ + (v^a • • • • 

The end point propagators are third kind polynomials, V n , while internal ones are 
second kind, U n . The determinant is 



J2 



h =1 j'i>j 2 =1 

( 75 ) 

+ E U Ji Ch ~ -72) ^2 C/2 - is) v j 3 + + ( v l v 2 ...V v ), 
3i>32>h = l 

since V n (x) is unity when A = 0. Det^ correctly vanishes when all the Va do due to 
the resulting zero mode. 

These expressions can be used to discuss a Borg-Levinson inverse theorem (e.g. 
Hald, [33]) and Dikii trace identities, but such further analysis and manipulations 
must be postponed. 
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